The imposition of Cauchy data to the Teukolsky equation III: The rotating case 
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We solve the problem of expressing the Weyl scalars ip that describe gravitational perturbations 
of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the 
Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express ip and dttp 
as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out 
of it) of a generic spacelike slice of the spacetime. These results provide the link between initial 
, data and ip to be evolved by the Teukolsky equation, and can be used to compute the gravitational 

■ radiation generated by two orbiting black holes in the close limit approximation. They can also be 

i—{ ' used to extract waveforms from spacetimes completely generated by numerical methods. 
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I. INTRODUCTION 



In Ref. |IJ the question was raised of how to impose initial data to the Teukolsky equation (that describe perturba- 
tions around a rotating black hole). We noted that the expressions of Chrzanowsky g] for the Weyl scalars V>4 and ipo 
_^ , in terms of metric perturbations were written as second order operators on the four-metric and appeared inconvenient 
^ ■ at the moment to use them for building up the initial values needed to start the integration of the Teukolsky equation. 

The work of reference Q| showed how to solve the problem for a nonrotating background, i.e. perturbations around 
a Schwarzschild hole by relating Weyl scalars vb, to the Moncrief waveforms 4>m, a- n alternative description of metric 
perturbations explicitly built up out of the three-metric 9^ and the extrinsic curvature Kij of the hypersurface t = 
constant. In Ref. || the ip - <j>M relations were successfully tested with a program for integration of the Teukolsky 
' equation. 

, It is not obvious how to extend the above techniques to the rotating case. Thus, in the present paper we turned to 
a more geometrical approach that lead us to the desired relations for rotating holes. In Sec. II we collect the results 
of the 3+1 decomposition reviewed in Ref. Q relevant for our derivation. This has the advantage that makes ip to be 
automatically independent of the shift, so our task is reduced to prove that terms depending on the first perturbative 
^Jy order lapse vanish. This is made in Sec. Ill, where we also build up dfip in terms only of Syand 2£y. This results 
allow to compare, given the initial data, evolution through integration of the full Einstein equations and Teukolsky 
equation (linearization around a Kerr hole), and test, for instance, the close limit approximation for orbiting holes. 
Notation: We use Ref. || conventions. An overbar on geometric quantities means that they are three-dimensional 

• quantities, i.e. defined on the t — constant hypersurfaces £4 (an exception to this rule is the complex conjugation 
of the vector m a , i.e. m Q ). (a, (3) and [a, (3] on indices a, (3 represent the usual symmetric and antisymmetric parts 
respectively. Greek letters indices run from to 3 while latin letters indices run from 1 to 3. Subindexes (0) and (1) 
mean pieces of exclusively zeroth and first order respectively. 



II. GEOMETRIC STRUCTURE AND GRAVITATION 



Following Ref. Q] we write the metric as 

ds 2 = -N 2 (e°) 2 +g ij e i e j , (2.1) 

with 6° = dt and l = dx % + N l dt, where iV 1 is the shift vector and N the lapse. 
The cobasis 8 a satisfies 

d0 a = --C^9 13 A I? 7 (2.2) 

with = —Cj Q = djN 1 and all other structure coefficients zero. Note that g~ij = gij and g 1 ^ = g % K 
The spacetime connection one-forms are defined by 
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where rjg 7 denotes the Christoffel symbol. These connection forms are written out explicitly in ||. In particular, 
lj 1 jk — r l jfe = T l jk , and the extrinsic curvature is given by 

K i j = -N^ j = ~N- l d g ij , (2.4) 

where we define the operator 

d = ^ - £ N , (2.5) 

with £n the Lie derivative on the hypersurface E t with respect to the vector N z . Note that do and di commute. 
The Riemann curvature tensor is given by Q 

R a f3 P a = dpU}0 a — d a uj a f3p + uj a \ p uj X i3 a — Lu a \ a uj X p p — uj a p\C x pa (2.6) 

For rewriting in the next section the Weyl scalars in terms of hypersurface quantities only, we relate the spacetime 
Riemann tensor components to the 3-dimensional Riemann and the extrinsic curvature tensors 

Rijkl = Rijkl + ZK^K^j, (2.7) 

R 0ijk = 2NVyK k]i , (2.8) 

RoiOj = N(doK l3 + NK ip K p 3 + ViVjN). (2.9) 
Another important relation in three dimensions is 

Rijki = ^fji[kRi\j + 2gj[iR~k]i + Rgi[Wk] r (2-10) 

The Ricci tensor R a p = R° acr p is given by 

Rij = Rij - N- 1 doK l3 + KK l3 - 2K lp K p J - N^ViVjN, (2.11) 
R oi =NW j {Kg i j -Kij), (2.12) 

Roa = N V 2 N - N 2 K pq K pq + NdoK. (2.13) 

In order to incorporate the source terms we consider the Einstein equations as R a p = T a p — ^g a pT. For instance, 
the "Energy constraint" is defined by 

G°o = ^(K mk K mk -K 2 -R)=T° . (2.14) 

Finally, from its definitions 

doRij = V fc (&r£.) - V.^offfc), (2.15) 

where 

&r£. = -2\7 {l (NK j} k ) + V k (NK i:j ). (2.16) 

Note that writing equations in terms of do instead of dt allowed us to get rid of the shift dependence. This is 
because do is orthogonal to the spacelike hypersurface S t . 

III. WEYL SCALARS FOR KERR PERTURBATIONS 

For the computation of gravitation radiation from astrophysical sources it is convenient to work with the Weyl 
scalar 

ipi = ~C a j3 1 sn a m l3 n < rfi 6 , 
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since it is directly related to the outgoing gravitational waves. For perturbations around a Kerr hole we have 



-tp4, — Rijkin l m? : n k m l 



4Ro j0 iri [0 m j] n [0 rnf ] . 



Eqs. ( |2.7| ) and ( |2.8| ) directly give us the two first terms in the above sum in terms of hypersurface geometrical 
objects (gij, Kij). In the last term we have to make use of Einstein equation ( 2.11 ) to eliminate d^K.^. If one now 
considers first order perturbations around a Kerr hole, one would have to consider in ipi two types of terms: terms 
that involve first order perturbative Ricmann tensors contracted with the background tetrads and terms that involve 
the Ricmann tensor of the background contracted with three background and one perturbative tetrads. It is not 
difficult to see that the latter terms vanish for the Kerr background. For the Kerr geometry the only non-vanishing 
Weyl scalar is ip2 = i? Q /3 7 5^ a m^n 7 m' 5 and one can quickly see that the above contributions, even with one of the 



tetrads being a perturbative one, still vanish. For instance, consider the term Rijkion^m 3 'n k 'rfv 



This term vanishes 

because it is contracted with two m vectors, and any contraction with a repeated tetrad vector of the Riemann tensor 
vanishes for the Kerr spacetime. Similar arguments apply to the other terms. 

Let us turn our attention to the ter ms t h at in volve the first order Riemann tensors contracted with the background 
tetrads. Taking a look at equations (p/7|)-(2.S) we see that if one considers first order perturbations, we will have 
expressions involving the first order extrinsic curvature, metric, and lapse. We do not want our final expression to 
depend on the perturbative lapse. It is easy to see that it actually does not depend on it. For i?oijfc we see that 
the lapse appears as an overall factor. So the expression evaluated for the perturbative lapse is proportional t o the 
expression evaluated in the background, which vanishes. For Roioj if we rewrite it using the Einstein equation (2.11) 
again the lapse appears as an overall factor and the same argument as for Roijk applies. As a separate check, we have 
verified the independence on the perturbative lapse and shift using computer algebra. 
The final result for the first order expansion of the Weyl scalar V>4 therefore is, 



1p 4 



Rijki +2K t [ k Kqj 



(i) 



n l m j n k m l 



47V, 



(0) 



K 



i[k,i\+ Tj[ k Ki]p 



+47V, 



(0) 



Rji -K 0P Kf + KK-ji - Tjt + -Tg 3l 



(i) 



n^m j ^n k m l 



(1) 



(3.1) 



where AT/ ) = G?ktrr)~ ^ s ^ ne zer °th order lapse, n l , m? are two of the null vectors of the (zeroth order) tetrad (see 
Ref. 0), latin indices run from 1 to 3, and the brackets are computed to only first order (zeroth order excluded). 
To obtain dtipA, the other relevant quantity in order to start the integration of the Teukolsky equation, we can 



operate with do on tp4 given by Eq. (3.1) to find 



d t ^ 4 = Nfadj (^ 4 ) - n l w?n k m: 



■J „ k m I 



doR 



ijkl 



(i) 



(3.2) 



+AN {0) n^n k m l d K j[k>l] + d T p j[k K l]p + f j[k d K l]p 



-ANl 0) n [0 m j] n [0 m l] 



(i) 

do Rji -2K^d K j)p - 2N i0) K jp K%K? 



-K 3l doK + KdoKjt - doTjt + - 9jl T - N [0) TK jt 



(i) 



where we made use of the equality 



>i,,iW = 2NK(. 



The derivatives appearing in Eq. (3^) can be obtained from Eq. ( 2.13 ) 



d K = N (0) K pq K pq - V N (Q) - N~]T, 



(0) 



(0) 



! 00, 



(3.3) 



fromEq. (2.14) 



do R= 2K^d K pq + m (0) K pq K p s K sq - 2Kd K - 29 T °, 



(3.4) 



and from Eqs. (2.10) and (2.4) 
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OoRvh = -4iV (0) {K i[k R t]j -K, [h - Rl]l -- R (K i[k9l]j - K^ m 

+2gi[kd Ri]j -2g j[k d R^ -g l [ k gi] j d R +2K i [ k d K l ] j - 2K j[k d K l] , 



(3.5) 



Note that in the last three equations we have taken explicitly the lapse to the zeroth perturbative order. This is 
so because in building up d t ipi explicitly all dependence on Nn\ cancels out. To prove this one can do the explicit 
calculation for the Kerr background using computer algebra. An alternative is to notice that <9o04 — <Ct"04 where 
t a is a vector that includes the background and first order perturbations of the lapse and shift. If one now expands 
out this expression one gets doV>4 — £t (0) V>4 (0) + Ct {0) ip4 {1) + £t (1) ^4(o) ■ Now, since ^4 (0) vanishes identically for all 
time, the only contribution one has is doip4 = £t (0) V^ij ■ Therefore the time derivative of ^4 does not depend on the 
perturbative lapse and shift, since neither £ t(0) (by construction) nor ip4 m (due to the proof we gave above), do. 

The other 
by Eqs. 



: her p ie ces n eeded to bu ild up d t ip4 only out of hypersurface data are d^Kij^oT^, and do Rij that are given 
( 2.1l| ), ( 2.16 ) and ( [2.15| ) respectively. As before, we have to consider the zeroth order lapse only, for instance 



d K i:j = iV (0) 



R i:j + KK i:i 



2K ip K p j - N^ViVjN (0) 



T, 



(3.6) 



(i) 



This completes our proof. A check of the relations (3T) and ( p.2[ ) can be made in the Schwarzschild background 
for close limit initial data where H at t=0 we have dttp = — "^fip- 



IV. DISCUSSION 



The issue of expressing ip explicitly in terms of hypersurface data only appears as of a purely technical character, 

but it is of great practical use. Especially when one thinks of the important role played by first order per turbat ion s 

as testbeds for comparison with full numerical integration of Einstein equations. Note that since Eqs. (O) and (gj) 

hold on any t — constant slice of the space time can not only be used to build up initial values for ip and dtip, but 

also at a later time to extract fully numerically generated wave forms. 

The above equations provides the desired link between initial data (consisting of g^ and Kij ) and the Weyl scalar 

-k - - 

-04. Geometrical objects like r^, Rij and Rij k i involve first and second derivatives of the metric. Since astrophysical 
i niti al data for Kerr perturbations are numerically generated H] this fact has to be taken into account. Expression 
(3.1) also includes a source term that allows to incorporate perturbations generated by particles or accretion disks 
around Kerr holes. 

If one chooses to work in the Teukolsky equation with %po = —C a ^&l a TnPl 1 m 1 ' , which gives a better representation 
of ingoing gravitational waves, a completely analogous procedure applies to connect it to hypersurface data upon 
replacement of the double contractions with the corresponding null vectors l a and mP instead of n a and m . 

Finally, we have been able to write ip^ and tpQ on the hypersurface £<, but we did not said why. In fact it is not 
warranted that one can do that with any object defined on the spacetime. Is this because they are first order gauge 
invariant objects? This shouldn't be enough since we checked that for tp 3 (and the same for tpi), we do not succeed 
in writing them in terms only of objects on the slice t = constant. The key point here seems to be that and ipo 
are also invariant under tetrad rotations and then directly connected to physical quantities, while ip3 ancl "01 are n °t- 
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APPENDIX A: ALTERNATIVE EQUATIONS 

We can put all this together to yield the following expression of the first order perturbation in ip4 in terms of 
perturbations in the 3-metric Sgij, perturbations in the extrinsic curvature SKij, and several quantites from the the 
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background (Kerr) geometry, the spatial metric ^g^y, the extrinsic curvature , the lapse and the shift 

. We have already argued that first order perturbations of the principal null vectors n M and rh^ will not contribute 
to Stp4 so we have 

Sip4 — 5A i jkin l m j n k m l + 2SB l j k n 3 m k [n°rh l — n l fh°] + SC t j[n m l n m j + n l m°n J fh° - n°m l n J m — n°fh?n l fh°] 
where 

SA ijkl = S ^R ljkl + [Kf5K ik + KjpSKfl - (k «-> I)] 

6B ijk = N^[D, 5K lk - \[D k S ^g ml + A S ^ g mk - D m S ^g lk ] i*)gW lm K \f - (k <- j)] 

K7« = AT (0)2 4%/ < 3 V ro + A^ 2 <L4 iym (3 V 0) ' m - [SBwNW + B^5^g lm N^\ m + A%5 ^g m N^ m 
+6A ja NW + 6Aa im N^ l N^ m + A(% jm N^ , k 6^g kl N^ m + A™ iljm N^ ? 6™ g *™ N (o) fc] 

and 

S < 3 >i$ w = \D h [ ^g^ m (D l S^g mj + D j 5^g ml - D m S^ 9jl )} - (fc <-> l) 

To calculate dtipA we use the above expression for 5ip4 and plug in dtS^gij and SdtKij for S^gij and SKij in the 
above, respectively. Where, dtS^gij and SdtKij can be obtained from Einstein's equations as follows: 

d t 8^g l3 = -W^SKv + N(°) k S^ gij , k + <W fc < 3 ty«\. fc + S^g^N^, 

J9 t ^. = I [A SWg mi + A <5 (3) 5mJ - An <y (s) 5«] V°\ fc 

+iV (°)[(5( 3 )i? y . - 2K^ k jSK ik - 25K k .jK^ lk + K^ i3 5K + K^SKy] 

+N^ k 5K l3 , k + 6K ik NW k j + 6K kj N^" ti + K<%[5& g^N™ k ] tj 
+K^ l3 [S^g kl N^ k i l + N^ l S^g lk K^^ k 

where SK = ^g^ lJ SK tJ + K^S ^ g J and 8K\ 3 = 8K jk ( 3 V°) fei + K^ jk S ^g H . 
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